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Abstract. We start with a Riemann-Hilbert Problems (RHP) 

WITH CANONICAL NORMALIZATION WHOSE SEWING FUNCTIONS DEPENDS 
ON SEVERAL ADDITIONAL VARIABLES. USING ZAKHAROV-ShABAT THEO- 
REM WE ARE ABLE TO CONSTRUCT A FAMILY OF ORDINARY DIFFERENTIAL 
OPERATORS FOR WHICH THE SOLUTION OF THE RHP IS A COMMON FUN- 
DAMENTAL ANALYTIC SOLUTION. THIS FAMILY OF OPERATORS OBVIOUSLY 
COMMUTE. Thus we are able to construct new classes OF INTE- 
GRABLE NONLINEAR EVOLUTION EQUATIONS. 

1. Introduction 

The development of the soliton theory revealed an important class of NLEE (non- 
linear evolution equations) that describe special types of wave-wave interactions 
[H [291 HI HI [231 [m [3D] which play important role in various fields in physics. 

A formal approach to the integrable equations started by Gel'fand and Dickey 
[21 [T7j and developed actively later on (see e.g. [T7] and the references therein) 
is well known. It allows one to construct the Lax representations for important 
classes of NLEE such as the dispersionless KP hierarchy but it disregards the 
spectral properties of the Lax operators. 



2000 Mathematics Subject Classification: 03.75.Lm, 42.65.Tg, 45.20.Jj 
Key words: Riemann-Hilbert Problems, commuting operators, integrable 3-wave interactions 



2 



y. S. Gerdjikov March 2, 2013 



The topic quickly attracted mathematicians from spectral theory, dynamical 
systems, Lie algebras, Hamiltonian dynamics, differential geometry, see [291 ESI 
[Ml ESI O [iHl H] and the numerous references therein. It attracted also a number 
of physicists because they found important applications of these NLEE in fluid 
mechanics, nonlinear optics, superconductivity, plasma physics etc. As a result 
many different approaches for investigating the soliton equations and construct- 
ing their Lax representations, soliton solutions, integrals of motion, Hamiltonian 
hierarchies etc. were developed, see [361 ESI El ESI ESI [3S1 EI]- Of course, it is 
not possible in a short paper to list all important references that cover the broad 
topics mentioned above. 

The inverse scattering method has been applied to many physically impor- 
tant multidimensional evolution equations including the N-wme equation, Davey- 
Stewartson, Kadomtsev-Petviashvilh etc. [Ml ED EH IMl ESI ESI EH] . They have 
been treated by nonlocal generalizations of the Riemann-Hilbert problem and by 
the 5-method. 

In the present paper we propose an alternative approach to the same class of 
equations using as a starting point the Riemann-Hilbert problem (RHP) [MIES 
[34 l l35 t[29ll37j : the importance of the canonical normalization of RHP was noticed 
in [101 [6] . Our aim is to show that this allows one to construct rings of commuting 
operators and in addition gives a tool to study their spectral properties. 

In Section 2 below we start with some preliminaries concerning the RHP. In 
Section 3 we use the solutions of the RHP to construct family of jets of order 
k, in Section 4 we list their simplest reductions. In the last two Sections we 
demonstrate how this construction can be used to solve NLEE in two and higher 
dimensional space-times. In Section 5 we use jets of order 1 to reproduce well 
known results about the 3-wave equations in two- and three-dimensional space- 
times. We also demonstrate the integr ability of N-wave type equations in higher 
dimensional space-times. In Section 6 we use jets of order 2 which allows us to 
construct new types of integrable A^-wave interactions whose interaction terms 
contain quadratic and cubic nonlinearities, as well as x-derivatives. These equa- 
tions also allow integrable extensions to three-dimensional space-time. The last 
Section contains discussion and conclusions. 

2. RHP with canonical normalization 

Let us formulate the RHP: 

(1) ^^{x,t,X) = Cix,t,X)G{x,t,X), A'^gM, lim ^+(x,t,A) = U, 
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where ^^{x,t,X) take values in the simple Lie group © with Lie algebra g. 
^~^{x,t,X) (resp. ^~{x,t,\)) is an analytic functions of A for ImA*' > (resp. 
for ImA*^ < 0). For simplicity we consider particular type of dependence of the 
sewing function G{x, t, A) on the auxiliary variables: 

(2) - X''[Js,G{x, t, A)] =0, i— - A'^ [K, G{x, t, A)] = 0. 

where A; > 1 is a fixed integer and Js are linearly independent elements of the 
Cartan subalgebra Js G f) C g. 

The canonical normalization of the RHP means that we can introduce the 
asymptotic expansion 



(3) C^(x,t,A) = expQ(f,t,A), Q{x,t, X) = J2Qk{x,t)X'' 



k=l 

Since ^^(x, t, A) are group elements then all Qk{x,t) G g. However, 

(4) Jsix,t,X) = C^{x,t,X)Jsi^ix,t,X), lC{x,t,X) =^^ix,t,X)Ki^ix,t,X), 

belong to the algebra g for any J and K from g. If in addition K also belongs to 
the Cartan subalgebra f), then 

(5) [Js{x,t,X),K:{x,t,X)] = 0. 

An important tool in our considerations plays the well known Zakharov- 
Shabat theorem [38l |39] formulated below 

Theorem 1. Let ^^{x,t,X) he solutions to the RHP ([Ip whose sewing func- 
tion depends on the auxiliary variables x and t via eq. Then ^^{x,t,X) are 
fundamental solutions of the following set of differential operators: 

Lst + Us{x, t, X)t{x, t, A) - X''[Js,C^{x, t, A)] = 0, 



(6) ± 

^i4r + ^(^' ^)^^ t , A) - A'^ [K, e± (f , t, A)] = 0. 



dxs 

Proof. The proof follows the lines of [6S\ I39j. We introduce the functions: 

gfix, t, A) = i^i^ix, t, A) + X'^^^ix, t, X)Jsi^{x, t, A), 



(7) ff^ 

g^{x, t, A) = i-lftix, t, A) + X'^tix, t, X)Ki^{x, t, A), 
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and using ([2]) prove that 

(8) gf{x,t,X)=g-{x,t,X), g-^ {x,t, X) = g' {x,t, X), 

which means that these functions are analytic functions of A in the whole complex 
A-plane. Next we find that: 

(9) lim gtix, t, A) = A^ J„ lim g+{x, t, A) = X^K. 

^ ' A^cxD A— J-oo 

and make use of Liouville theorem to get 



gt{x,t,X)=g-(x,t,X) = X^J,-^\J,,i(x,t)X^~\ 

1=1 

k 

g+{x, t, A) = g~{x, t, A) = X^K - ^ Vl(x, t)A^-'. 



1=1 

We shall see below that the coefficients Us-i{x,t) and Vi{x,t) can be expressed in 
terms of the asymptotic coefficients Qs in eq. ([3|). □ 

Lemma 1. The set of operators Lg and M commute, i.e. the following set 
of equations hold: 

i^ - + [Us{x, t, A) - X'^Js, Uj{x, t, A) - X'^Jj] = 0, 

, , OXj OXg 

^ ^ du dv 

where 

k k 

(12) Us{S, t,X) = J2 i)A^-', V{x, t,X) = Yl ^(^' ^)^''"'- 

1=1 1=1 

Proof. The set of the operators Lg and M Q have a common FAS, i.e. 
they all must commute. The eqs. (jlip are an immediate consequence of (l6|). □ 

3. Jets of order k 

In what follows we will consider the jets of order k of J'(x,X) and /C(x,A), see 
([5]). We introduce them by: 

Js{x,t,X) = (x^C^ix,t,X)Jii^{x,t,X)) = A^J, -C/,(f,t,A), 

(13) ; ^ 

}C{x, t, A) = (A^e^(f , t, X)K^^{x, t, A) j = X'^K - V{x, t, A). 
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The subscript + used above means that we insert the asymptotic expansions of 
and their inverse ([3]) and cut off the terms with negative powers of A. 
Obviously Us{x) G g can be expressed in terms of Qs{x)- In doing this we 

take into account (El) and obtain [T8] 



oo ^ oo ^ 

(14) J,(f,t,A) = J, + ^^ad^J„ }C{x,t,\) = K + Y-ad''QK, 

k= 

and therefore for Us-i we get: 



k=l k=l 



Us-i{x,t) = -adg^Js, Us-2{x,t) = -adQjJs - ^adg^J^ 

1 1 o 

Us;3{x,t) = -adq^Js - - (adg^adg^ + adg^adgj Js - ^adg^^ 



(15) ■ ^ 

Us-k{x,t) = -adg^^ Js - - ^ adg^adQj^J^ 



2 

s+p=k 



~ ^ X] adQ^adQ^adg, -^ad^^^, 

s+p+r=k 

and similar expressions for Vi{x,t) with Jg replaced by K. 



4. Reductions of polynomial bundles 

An important tool to construct new integrable NLEE is based on Mikhailov's 
group of reductions [22]. Below we will use mainly Z2 and Zjv with N > 2 
reduction groups. The basic Z2-examples are as follows: 

a) AS,'^'\x,t,€\*)A = C{x,t,X), AQ\x,t,€\*)A = -Q{x,t,X), 

(16) b) B^+'*{x,t,eX*)B = r{x,t,X), BQ*{x,t,eX*)B = Q{x,t, X), 
c) CC+^^{x,t,-X)C = r{x,t,X), CQ^x,t,-X)C = -Q{x,t,X), 

where = 1 and A, B and C are elements of the group S5 such that = B^ = 
= U. As for the Z^v-reductions we may have: 

(17) DC^{x,t,u}X)D = C^{x,t,X), DQ{x,t,ojX)D = Q{x,t,X), 

where = 1 and = U. 

These relations allow us to introduce algebraic relations between the matrix 
elements of Q{x,t,X) which will be automatically compatible with the NLEE. 
The classes of inequivalent reductions of the iV-wave equations related to the 
low-rank simple Lie algebras are given in [51 [9l [TTl [12l [131 IIS] • 
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5. On A^-wave equations {k = 1) in 2 and more dimensions 

The integrability of the N-wave equations has been well known for several decades 
now, f25l [Ml [Ml [SSI El 1211 [22]. Their Lax representation involves two Lax 
operators linear in A which are particular case of ([6|) with k = 1: 



(18) 



dx 



+ [J, Qix, t)]e^(x, t, A) - A[J, e^(x, t, A)] = 0, 
+ [K, Q{x, {x, t, A) - X[K, e±(x , t, A)] = 0. 



The corresponding equations take the form: 



(19) 



dt 



ox 



[[J,Q],[K,Q{x,m=0 



In fact the construction of the FAS for the operator L (USD [351 IB^l [5U1 [25] was 
the important step forward, that demonstrated the importance of the RHP for 
solving integrable equations. 

The most important and nontrivial example of such NLEE is the 3-wave 
equations in two-dimensional space-time [361 129j- The most important and non- 
trivial case corresponds to g ~ sl{3) 



(20) 



Q{x,t) 



ui Us 

-Vi U2 
-V3 -V2 



J = diag(ai, 02,03), 
K = diag (61, 62,^3), 



with tr J = tr X = and oi > 02 > 03. We also impose the reduction (116k ) with 
A = diag (1, ei, £2) where e\ = e\ = 1. Then the 3-wave equations take the form: 



-I- Keie2ii2ii3 = 0, 
-|- Keiii*U3 = 0, 
-|- Ke2ulu2 = 0, 



dui 


ai - 


a2 dui 


~dt ~ 


61- 


62 dx 


du2 


a2 - 


as du2 


~dt ~ 


62- 


63 dx 


du3 


ai - 


as du3 


at 


bi- 


63 dx 


= ai(&2 


-h) 


- a2{b 



(21) 

where 
(22) 

Depending on the choice of the reduction and on interrelations between the group 
velocities the 3-wave interactions may describe qualitatively different processes: 
soliton decay and explosive soliton instability [Ml [29] . 
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In the case of 3-dimensional space-time we consider Q of the form (|20p . but 
now let Uj and Vj be functions of xi = x, X2 = y and t. Let also Ji = J and 
J2 = I = diag (ci, C2.C3). Now the corresponding solution of the RHP ^^(x, y, t, A) 
will be FAS not only of L and M above, but also of 

(23) + [/, Q(x, t)]C+(x , t, A) - A[/, ^+(5^, t, A)] = 0, 

and all these three operators will mutually commute, i.e. along with [L, M] = 
we will have also [L, P] = and [P, M] =0. As a result Q{x, y, t) will satisfy two 
more NLEE of the form (j24p . Obviously it will satisfy also 



dui ai-a2dui ai - 02 Sui ^ 

~ T ^ — H 5~ + ('^1 + '^2)eie2U2U3 = 0, 

OT oi — 62 ax ci — C2 ay 

ot 61 — 03 ax ci — C3 ay 
5m3 02-039^3 02-035^3 ^ ^ 

at 62 - 03 ax C2 - C3 dy 
which is linear combination of the three equations mentioned above. Here ki = k 
(see eq. (j22|) and 

(25) K2 = ai(c2 - C3) - 02(01 - C3) + 03(01 - C2). 

These three wave equations are related to the real forms of the algebra sl{3) 
which has rank 2. Therefore, trying to add more auxiliary variables to the solution 
of the RHP will not be effective since only two elements of all XgJs will be linearly 
independent. 

For N-wave equations related to Lie algebras q of higher rank r we can add 
up to r auxiliary variables. The corresponding PDF takes the form: 

(26) - E(ad7>dj)^ - i ^^j} U Q], [Js,Qix, t)]] = 

s=l ^ s=l 

where Q is an n x n off-diagonal matrix depending on r + 1 variables. We remind 
that if J = diag (oi, . . . , a„) then 

(adjQ)jfe = ([J, Q])jk = (oj - ak)Qjk, iadj^Q)jk = Qjk, 

Uj — Ofc 

and similarly for the other Jg. The coefficient r multiplying the i-derivative can 
be removed by rescaling of t. 

Again we can use additional reductions of the type (I16p . More details about 
these equations will be given elsewhere. 
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6. New A^-wave equations (k = 2) in 2 and more dimensions 

Here we shall give examples of new types of A^-wave equations. Let us choose 
again g = s/(3). The general form of the potentials is given by 

/ -ui ^3 \ / qii wi W3 

(27) Ql{x,t) = -Vi U2 \ , Q2{x,t) = -Zi q22 W2 

\ -V3 -V2 / \ -Z3 -Z2 qsS 

We also fix up A: = 2. Then the Lax pair becomes 

Le± = + U{x, t, A)e^(x, t, A) - AV, t, A)] = 0, 



(28) 

Mt = + V{x, t, A)^±(x, t, A) - X^]K, ^^(x, t, A)] = 0, 
where using eq. (fT5|) 

U = U2 + XUi = (^[J,Q2{x)]-^[[J,Qi],Qiix)]^ +A[J,gi], 

V = V2 + XVi = (^[K, Q2ix)] - ^[[K, Qi], Qi(x)]^ + X[K, Qi]. 

Note, that this Lax pair is independent of the diagonal elements of Q2- 

If we retain the generic potentials ()27p the Lax pair above will provide us 
with a set of 6 new complicated equations for the 6 independent functions uj and 
Wj. To make the things more simple we impose a Z2-reduction of the form (jl6b .) 
with A = diag (1, e, 1), = 1. Thus Qi and Q2 get reduced into: 

-ui \ / W3 

(29) Qi = \ eu\ U2 ] , Q2 = 

eu* / \ tyj^ 

and J and K are as in ()20p . Now L and M involve only 3 independent functions. 
Skipping the details we get a new type of integrable 3-wave equations: 

i[ai - a2j^r- - t{bi - h2)-z h eKUstis + e- ^ui\u2\ = 0, 

ot ox (ai — as) 



(30) 



i[a2 - 03)—- - i(62 - 03)^:^ + eKu^u^ - e- -\ui\ U2 = 0, 

ot Ox (ai — 03) 

5ii3 (9^3 m d{uiU2) 

^^"^ - - ^^'^ - 5^ - ai-a3 dx 

(ai-a2. |2 ,02 - 03 2\ , n |2 I i2\ n 

Vai - 03 ai - 03 y 
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where the interaction constant k is given by (j22p and: 

2o2 — ai — 03 

(31) U3 = W3 + — ^UiU2. 

2{ai - 03) 

The diagonal terms in the Lax representation are A-independent. Two of 
them read: 

H^l - «2)^: HOl - 02)^:^ eK(Mi'U2M3 - U1U2U3) = 0, 

/Qr,\ at ox 

^ ) ■^\ i2 ■^\ 12 

Z(a2 - «3)— HO2 - eK(UiU2U3 - U1U2U3) = 0, 



These relations are satisfied identically as a consequence of the NLEE (j30p . The 
third one also vanishes since tr [L, M] = 0. 

Let us now consider the case when the sewing function G of the RHP depends 
on 3 variables: t, xi = x and X2 = y with Ji = J and J2 = / = diag (ci, 02,03). 
For k = 2 we, obtain a set of three ordinary differential operators: L, M (I28p and 



= + y, t, A)e^(x, y, t, A) - A^]/, ^^(2:, y, A)] = 0, 
(33) W = W2 + XWi 

[/, Q2(x, y, t)] - ^[[/, Qi] , Qi(x, y, t)]^ + A[/, Qi (x, y, t)] , 

commuting identically with respect to A. It is obvious that [i,-P] = if 

■( sdui , * , '^2(01 - 02) I ,2 n 

z(ai - 02)^7- - t{ci - C2)^ h eK2U2U3 + e— n2 = 0, 

dt dy (ai - 03) 

V x5?^2 \9U2 ^ ^ ^2(02- as) I |2 r, 

i[a2 - 0.3)^- - i[c2 - C3)— h eK2UiU3 - e— ^\ui\ U2 = 0, 

.g^x ot dy (ai - 03) 

5m3 9w3 iK2 d{uiU2) 

nai - - nci - C3)^ — 

at oy ai — 03 oy 

(CL\ — 0,2 |2 ^2 — 13 I i2 \ /I i2 I i2\ 
H |n2| nin2 + eK2'U3(|'Ui| - P2I ) = 0, 
oi - as ai - a3 / 

where K2 is given by eq. ()25p . It is not difficult to write down the third new 
3-wave equation which is a consequence of the commutation [M, P] = 0. 

Since the three operators L, M, and N mutually commute, ui, U2 and n3 
as functions of x, y and t should satisfy simultaneously the three NLEE of the 
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type (j30p . Therefore they should satisfy also any NLEE which is obtained as, say 
linear combination of the above: 

2i— - - t{vn) ■ Vjlii + e{Ki + K2) h 7 r = 0, 

ot ^ ' \ai - a2 (ai - as) J 



(35) 



2z— - - z(w(2) • V)'U2 + e(Ki + K2) — = 0, 

ot \ai - as (ai - as) / 

2z— - z(v(s) • V)u3 - z— — + (Iml - \u2\ )ns 

ot ^ ' (ai — as)^ ai — 03 

+ 7 yT ((Ol - 02)^1 + (02 - as) M2 ) U1U2 = 0. 

(ai - as)^ 

Here V = {dx,dy)^, the characteristic velocities vq), j = 1,2,3 and k are two- 
component vectors given by: 

1 f h-b2\ . 1 

^ ^ ai - a2 V ci - C2 y ^ ' 02 - as V C2 - cs 

(36) 



-'(3) 



and Ki = K, see eq. (p2|) . 



61 


-62 


Cl 


- C2 


&1 


-63 


Cl 


- C3 



't2 



7. Discussion and conclusions 

We have proposed a method for constructing families of commuting operators. 
Applied to jets of order 1 with g ~ sl{3) this method reproduces the well known 
results for the 3- wave equations in two- and three-dimensional space-times. It is 
shown that A^-wave equations related to Lie algebras of rank r allow integrable 
extensions to r + 1-dimensional space-times. Below we briefly discuss some open 
problems and generalizations. 

Using jets of order 2 gives us the simplest nontrivial examples for new types of 
integrable 3-wave equation whose interaction terms contain quadratic and cubic 
nonlinearities, as well as x-derivatives. These equations also allow integrable 
extensions to three-dimensional space-time. 

It is not difficult to obtain many other new integrable 3- and A^-wave equa- 
tions. Indeed, one can choose: i) higher rank simple Lie algebras; ii) different 
types of grading; iii) different power k of the polynomials U {x, t, A) and V{x, t, A) 
and iv) different reductions of U and V. 

These new NLEE must be Hamiltonian. It is natural to view the jets U{x, t, A) 
as elements of more complicated co-adjoint orbits of the relevant Kac-Moody 
algebra, generated by the chosen grading of f, see [Ml [STj [32]. 
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By construction, the method allows treating multi-dimensional NLEE. In the 
examples above we used the algebra sl{3) and demonstrated integrable 3- wave 
equations in 2 + 1-dimensional space-time. If we want to study new types of 
integrable N-wscve models in r -|- 1 space-time dimensions we have to consider Lie 
algebras of rank r and accordingly larger values for A'^. 

The method allows one also to apply Zakharov-Shabat dressing method [351 
[39 | I27 1 [20] for constructing their explicit (A^-soliton) solutions. Instead of solving 
the inverse scattering problem for L we would rather deal with a Riemann-Hilbert 
problem with canonical normalization. For polynomials of order k the contour 
on which the RHP is defined consists of k straight lines Ik- argA = ni/k passing 
through the origin. Of course, it may necessary to use dressing factors with more 
specific A-dependence. 

This approach can be used also to analyze the NLEE derived by Gel'fand- 
Dickey approach [U [I7j . It would provide the possibility to systematically con- 
struct the spectral decompositions that linearize the relevant NLEE [71 |T6] . Still 
more challenging is to study the soliton interactions of the new N-wave equations. 
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